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We investigate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction arising from the hyperfine coupling 
between localized nuclear spins and conduction electrons in interacting carbon nanotubes. Using the Lut- 
tinger liquid formalism, we show that the RKKY interaction is sublattice dependent, consistent with the spin 
susceptibility calculation in noninteracting carbon nanotubes, and it leads to an antifen'omagnetic nuclear spin 
helix in finite-size systems. The transition temperature reaches up to tens of mK, due to a strong boost by a 
positive feedback through the Overhauser field from ordered nuclear spins. Similar to GaAs nanowires, the 
formation of the helical nuclear spin order gaps out half of the conduction electrons, and is therefore observable 
as a reduction of conductance by a factor of 2 in a transport experiment. The nuclear spin helix leads to a density 
wave combining spin and charge degrees of freedom in the electron subsystem, resulting in synthetic spin-orhit 
interaction, which induces non-trivial topological phases. As a result, topological superconductivity with Majo- 
rana fermion bound states can be realized in the system in the presence of proximity-induced superconductivity 
without the need of fine tuning the chemical potential. We present the phase diagram as a function of system 
parameters, including the pairing gaps, the gap due to the nuclear spin helix, and the Zeeman field perpendicular 
to the helical plane. 

PACS numbers: 71.10.Pm, 74.20.-z, 75.70.Tj, 75.75.-c 


I. INTRODUCTION 

The search for topological superconductivity and exotic 
quasiparticles supported by it, such as Majorana fermions 
(MFs), remains an ongoing challenge. MFs, being their own 
antiparticles, exhibit non-Abelian statistics and are promis¬ 
ing candidates for realization of quantum computationJ*^ De¬ 
spite intensive experimental^ and theoretical^ril efforts, the 
observation of MFs still remains inconclusive. It is there¬ 
fore important to propose experimentally achievable devices 
for the realization of such particles and even more exotic 
parafermionsji^— In this respect, carbon nanotubes (CNTs) 
seem promising due to advantages such as the availability 
of high-quality samples, high tunability, and strong electron- 
electron interactions due to the spatial confinementrSiS^ 
which is crucial for fractional statistics. 

In a nanowire with free carriers (electrons) and localized 
spins, such as spins of atomic nuclei, these two subsystems 
are coupled by the hyperfine interactioni2ii2^ With parame¬ 
ters typical for semiconducting nanowires, this interaction is 
weak on the scale of the electronic Fermi energy. It can then 
be recast as the Ruderman-Kittel-Kasuya-Yosida (RKKY) ex¬ 
change interaction^Si^ the electron mediated pairwise inter¬ 
action between localized spins (see also Ref. for systems 
beyond the RKKY picture). The strength of this pairwise in¬ 
teraction is given by the many-body state of the electron sub¬ 
system, and reflects its properties. For example, the RKKY 
coupling as a function of distance is modulated at the electron 
Fermi wavelengths^ while the spin-orbit interaction of elec¬ 
trons results in spin anisotropies^^*^ or suppression^ of the 
RKKY coupling. In low dimensional systems, the effect of 
electron-electron interactions becomes striking.i^^^ Namely, 
the RKKY interaction is strongly enhanced, formally seen 
as the renormalization of the exponent describing the power- 


law decay of correlators calculated in the Luttinger liquid for- 
malismi^ The stronger the electron-electron interactions, the 
more pronounced is this enhancement around the electronic 
Fermi momentum, which leads to a sharp resonant peak in the 
RKKY coupling in one-dimensional systems. At low enough 
temperature, this peak results in the ordering of the localized 
spins into a helix, which corresponds to an effective rotating 
magnetic held (Overhauser held) seen by the electrons. This 
macroscopic held changes the electronic state by opening a 
partial gap at the Fermi energy. While this further boosts the 
RKKY coupling strength, it is also interesting on its own. Be¬ 
cause a rotating held has a definite helicity, the partial gap 
opens in a spin selective way^Sill and the electronic subsys¬ 
tem also becomes helical. It has been theoretically suggested 
to exploit such helical Overhauser fields for, e.g., dynami¬ 
cal nuclear spin polarization^^ stabilization of fractionalized 
fermions^^ or production of tune-free topological matter4i^ 
Signatures of such a partial gap opening have been observed in 
GaAs quantum wires in transport experiments at sub-K tem- 
perature4i Subsequent density matrix renormalization-group 
analysis^ also supported the formation of the RKKY-induced 
magnetic order discussed in Refs. ljTlI^ 

In this paper we revisit the above picture considering metal¬ 
lic CNTs enriched by the atomic isotope with nuclear 
spin While Refs. also considered CNTs, the 

presence of sublattices was omitted. On the other hand, the 
results for the spin susceptibility calculated in the noninter¬ 
acting limit^i^T^ suggest that the RKKY interaction is lo¬ 
cally (between nearest neighbors) antiferromagnetic, unlike 
in, e.g., GaAs where it is ferromagnetic. The question there¬ 
fore arises whether in the presence of strong electron-electron 
interactions, the RKKY interaction retains its locally antifer¬ 
romagnetic character, and whether a macroscopic Overhauser 
field can still arise, which is necessary to push the transition 
temperature of the nuclear order to experimentally achievable 
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values, and to offer enriched CNT as a self-tuned topo¬ 
logical matter platform. 


To this end, we derive here the RKKY interaction taking 
into account sublattices explicitly in a CNT with interacting 
electrons. We find that the RKKY interaction is sublattice 
dependent, consistent with Refs. [34ll53l - lY^ and it leads to a 
locally antiferromagnetic nuclear spin helixi^ However, de¬ 
spite a lack of macroscopic spin polarization, the helix tran¬ 
sition temperature is still strongly enhanced, and reaches sev¬ 
eral tens of mK. We also confirm that the nuclear spin helix 
combining charge and spin degrees of freedom generates syn¬ 
thetic spin-orbit interaction for electrons, suitable to induce 
non-trivial topology^iiSiii supporting without involv¬ 

ing intrinsic spin-orbit interactions that happen to be still weak 
in CNTsi^— Therefore, we suggest to pursue experimentally 
the possibility to establish the RKKY induced nuclear spin 
order at low temperatures in CNTs highly enriched by 


This paper is organized as follows. In Sec.|II]we obtain the 
RKKY interaction within the Luttinger liquid formalism: in 
Sec. Ill Al we first establish the RKKY Hamiltonian in terms 
of the spin susceptibility of the conduction electrons, which 
are described as Luttinger liquid in Sec. Ill B1 the bosonization 
of the electron spin operators, which enter the spin suscepti¬ 
bility, are discussed in Sec. Ill Cl finally, using the results of 
Sec. Ill Ami cl we obtain the RKKY interaction in Sec. Ill PI 
The resulting Hamiltonian then represents a spin model for 
the nuclear spins, allowing us to investigate the nuclear spin 
order in Sec. [Till we take the ansatz for antiferromagnetic nu¬ 
clear spin helix in Sec. IIII Al and compute the magnon (spin 
wave) spectrum; in Sec. IIII Bl the transition temperature of 
the nuclear spin helix (without the feedback) is estimated. In 
Sec. IIVI we examine the feedback effect due to the nuclear 
Overhauser field: in Sec. IIV Al we show that the intervalley 
back scattering terms in the spin susceptibility are suppressed 
by the feedback; in Sec. IIV Bl we analyze the renormalized 
Overhauser field; in Sec. IIV Cl we estimate the transition tem¬ 
perature in the presence of the feedback, which is enhanced 
by more than four orders of magnitude. In Sec.|V]we investi¬ 
gate how the proximity-induced superconductivity affects the 
RKKY interaction: in Sec. IV Al we compute the spin suscepti¬ 
bility in the presence of the pairing gap; in Sec. IV Bl we show 
that the reduced transition temperature may still be within ex¬ 
perimentally accessible regimes. In Sec. [VI] we focus on the 
topological properties: in Sec. IVI Al a refermionized Hamil¬ 
tonian is established, which allows us to find MF solutions 
straightforwardly; the topological phase diagram is presented 
in Sec. IVI Bl Finally, we give a discussion on the nuclear spin 
helix and MFs in CNTs in Sec. I VIII The details of the calcu¬ 
lations on the spin susceptibility in the presence of the pairing 
gap and solving the Schrodinger equation for MF solutions 
are given in Appendix lAl and iBl respectively. 


II. R KK Y INTERACTION 
A. Hyperflne and RKKY Hamiltonians 

Nuclear spins of atoms embedded within CNTs couple 
to conduction electrons via the hyperfine exchange interac¬ 
tion. We consider a single-wall armchair-edged nanotube^ 
with the Hamiltonian 

= + Uu- ( 1 ) 

Here 'Heh discussed in Sec. IIIBI describes the interacting 
conduction electrons, and 'Hhf is the hyperfine coupling be¬ 
tween the conduction electrons and localized nuclear spins. 
The dipolar interaction between the localized spins is much 
smaller than these two terms and hence neglected.— (for sys¬ 
tems where this is not the case, the combined effect of di¬ 
rect and RKKY interactions may lead to, e.g., a canted spin 
state—) 

Assuming the electrons are in the lowest transverse mode 
due to a large transverse level spacing of the order of eV, we 
obtain an effective one-dimensional hyperfine interaction, 

= ( 2 ) 

where a = A, B denotes the sublattice index, j = 1,... ,N 
is the site index of cross sections along the tube axis, and 
Aq is the hyperfine coupling constant. There is a discrep¬ 
ancy between the measured hyperfine coupling constant and 
the theoretical prediction. The observed value in an isotopi- 
cally enriched (^ 99% ^^C) nanotube quantum dots^i*^ was 
two orders larger than the theoretical calculation employing 
a noninteracting system calculation—^ Whereas the mea¬ 
sured value was extracted through theories developed for other 
materials without valley degrees of freedom,— such as GaAs, 
and needs to be further confirmed)^ we take Aq = 6.0 /reV, 
which is in the order between the observed and theoretical 
values, for the purpose of estimation. We also note that 
subsequent measurements in CNTs with natural abundance 
(^ 1% i3c)^2i2S corroborate the hyperfine coupling constant 
reported in Refs.jsHH^. 

We split the nanotube into small cylinders of height a, the 
length of the atomic scale, and get N±_ = irRani as the num¬ 
ber of the atoms on each sublattice in such a cylinder, with 
71/ the atomic area density of a graphene sheet. We group 
the nuclear spins within one cylinder into an effective com¬ 
posite spin, iairj) = Sj/l which we refer to as a 

spin of a cross section. Because N± ^ 1, the effective spins 
are large, with maximal magnitude N±I, and thus can be 
treated semiclassically. We choose a to be the lattice constant 
of the CNT (the carbon-carbon bond length is a/s/3), which 
for (6,6) CNTs gives Nj_ = 12 and the radius R ^ 4.1 A. 
The choice of (6^6) CNTs is partially motivated by the ex¬ 
periment in Ref.j^ which reports that defect-free CNTs with 
definite chiral index have been made possible. We also de¬ 
note the nuclear spin of ^^C as / = 2. Further, the effective 
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one-dimensional electron spin density operator is 

Sa(r) = — - n) CTj, (3) 

n/ 

i 

where r denotes the coordinate along the tube, the posi¬ 
tion operator of the ith electron, cr^ is a vector with compo¬ 
nents formed by the Pauli matrices in spin space of the ith 
electron, and 4’^,a is the transverse part of the electron wave 
function (assumed to be the same for all electrons). We will 
assume 4^,a spreads uniformly over the circumference, so 
that \4^.,a'^ = 

Since Aq ^ ep, we use the Schrieffer-Wolff transforma¬ 
tion to integrate out the electron degrees of freedom ,^^'^^'^^ 
which results in an effective RKKY interaction between two 
localized spins. 


'Hrkky = ^ 

where p, i/ = x,y,z are coordinates in spin space and the 
effective RKKY exchange coupling, 

-rj) = -r,). (5) 


electrons are confined into the lowest transverse mode due to 
the spatial confinement, the tight-binding model results in^i 


(^i7.a(9) 4.7 .<t(9)) 

( 0 -'yhvpq 

^ —"fhvpq 0 


,a{q') 


, (8) 


where cj, ^ „{q) is the creation operator with the sublattice in¬ 
dex a = A, B {a = ± 1 ), valley index 7 = ±, spin a =t, 
the z component of the momentum q = q^ is measured from 
the Dirac point K.^, and vp = is the Fermi velocity. 

Eq. ([ 8 ]l can be diagonalized by symmetric (S = -f) and anti¬ 
symmetric {S = —) combinations. 


4s,l,<r{q) — \cA,~f,a{q) + 5cB,~i,cr{q)\ , (9) 


corresponding to the eigenvalues, Es-y = —Sjhvpq. There¬ 
fore, the energy spectrum of Eq. ® exhibits linear dispersions 
close to the Dirac points, leading to two copies of Luttinger 
liquid spectrum located at kz = Eky (see Eig.[Tli. 

To proceed, we describe the system in terms of the right 
(R = -fl) and left (L = —1) moving particles, 4e,'y,<j{q)^ 
where for £ = R and L, we have 


The static spin susceptibility is defined as 


Xa'isin-rj) =/ die {[S^inR), S^rj ,0)]) , 


(6) 


with an infinitesimal positive rj and (...) being the average cor¬ 
responding to the one-dimensional effective electron Hamil¬ 
tonian Tfei- Eor the continuum description, we will replace 
S^{ri,t)/a with the operators (r, f). 


B. Electron Hamiltonian and bosonization 


4R,iAq) = '^saA<i)\s=-~i 1 (lOa) 

4L,j,a{q) = 4s,-f,a{q)\g—^ 1 ( 10 b) 

respectively. Combining Eqs. (|9l) and ( fTOl i. we find the rela¬ 
tion between the original electron operators and the right/left 
movers. 


t-a, 7 ,cr (^) 


1 
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4i,i,<y{q)\e=-y + oi4i,i,<y{q)\e=-y ■ 

( 11 ) 


One may bosonize 4i,-y,rT in terms of the bosonic fields 9^^ 
and 4'yrT In real space, we have 


In this section, we discuss the one-dimensional effective 
electron Hamiltonian and its bosonized form. We start with 
the Hamiltonian of an interacting electronic system, 

Ua =^-0+ Tfint, (7) 

where "Ho int describe the kinetic energy and interaction 
terms, respectively. 

The Hamiltonian "Hq is defined by a tight-binding model 
of a carbon lattice, including the nearest-neighbor hopping 
terms with the hopping parameter t. We neglect the longer- 
range hopping, nanotube curvature, and spin-orbit interac¬ 
tions, which results in a Hamiltonian conserving the total spin. 
The spin susceptibility can then be written as — rj) = 

~ Eourier transform "Ho expand it 

around the Dirac points, K.y = 'ykyZ + with ky = ^ {z 

and t being the unit vectors along the tube axis and circumfer¬ 
ence, respectively)^ With the assumption that the conduction 


4t,iAA 


^i(ekF+7ky)r^-i\e4>^^(r)-e^^(r)] ^^ 2 ) 

y/2Tra 


where kp is the Eermi wave number, the lattice constant a 
sets the smallest length scale of the system, is the Klein 

factor removing a {£, 7, a) particle from the system, and the 
bosonic fields satisfy the following relations^^ 

TT 

= i-sign{r 2 - ri)S^YSaa', (13) 

V 4 ^a{r) = -TT [pR,^,a{r) + PL,j,a{r)] , (14) 
V9^„{r) = TT [pR.j^a{r) - PL,'f,<j{r)] , (15) 

with the real-space density operator pi,-yAA ~ 
'4’i-y One can see that the field S/O^ctA/t: is 

canonically conjugate to 


4'icr {fl ), 


TT 


iSA - ri)Sjj4 aa' • 


(16) 
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Index 

Degree of freedom 

Possible values 

a 

sublattice 

+ 1 ),B(= 1 ) 

7 

valley 

-f, - 

(j 

spin 

t (=+!),;(= - 1 ) 

5 

symmetric/antisymmetric 
combination of a = A, B 

-f, - 

£ 

right/left mover 

R(=-bl),L(=-l) 

V 

charge/spin sector 

c(=-bl),s(= - 1 ) 

p 

symmetric/antisymmetric 
combination of 7 = ± 

S(=-bl), A(= - 1 ) 


TABLE I. The indices defined in Sec. mi the corresponding degrees 
of freedom, and the possible values of the indices. 


Including the electron-electron interaction, the electron 
Hamiltonian Tfei can then be bosonized^^ 

'Hei = ^ f — —[V0i,p(r)]^ 

+ J^ [V<ii',.p(r)]^} , (17) 

where i/€{c=-|-,s = —} refers to the charge/spin sectors, 
and P G {S =+, A = —} the symmetric/antisymmetric 
combination of the bosonic fields between the 7 = ± valleys, 
namely, 

0vP = 2 [^+4 + ^^+4 + P (^-4 + ^^- 4 )] I ( 1 ^) 
\ + v(j)+.x P P + v(t>-^i)]. (19) 

The velocities for the {i',P) channels are itj^p = vp/K^p. 
The indices defined in this section are summarized in Table U 
for reference. The noninteracting case corresponds to the Lut- 
tinger liquid parameters K^s = Kca = = Ksa = 

and the repulsive electron-electron interaction leads to KcS < 
1. The parameter KcS depends on the radius of CNTs through 

the relation, KcS = [l + (8e^)/(7r/iz;p) In (i?s/i?)] 
where e is the electron charge and Rg ~ 1000 A is the 
screening length.— Therefore, for CNTs with smaller radius, 
the electron-electron interaction has stronger effects due to the 
stronger spatial confinement, as expected. However, this ra¬ 
dius dependence is relatively weak because of its logarithmic 
form. For R = 4.1-100 A, KcS ~ 0.16-0.24. In this work, 
we take K^s ~ 0.2 and K^a « Kgg « K^a ~ 

With Eqs. (fT^ . ( fTSb , and ( fT9l l, we can write the single¬ 
particle spin operator (r) in terms of the bosonic fields to 
compute the correlation functions in Eq. (|6]l. Since the elec¬ 
tron Hamiltonian (fTTl i is a free bosonic system, the correlation 
functions can be computed straightforwardly within the Lut- 
tinger liquid formalism^ 


C. Spin operator in terms of the bosonic fields 

To examine the sublattice dependence, we first write the 
spin operator in terms of the original electron operators with 


the explicit sublattice index a. 


SaiTj) = (20) 


cr,fT' 7,7' 


which, according to Eq. (HHl, can be written as Sl^{rj) = 
where 


arises from the forward scattering (g ~ 0 or g ^ 2fc„) and 


,( 21 ) 


^,7 




,( 22 ) 


corresponds to the back scattering [g ~ 2 kp or g ^ 2(fc„ ± 
kp)]. Here are the Pauli matrices in spin space, i = —£, 
7 = — 7 , and the inverse Eourier transform of ^( 9 ) is 
given by 


Y (23) 

with defined in Eq. ( fTOl i. 

Since we consider the temperature T much lower than the 
Eermi energy ep, the states below ep are filled, allowing us 
to keep only the scattering processes that take place on the 
Eermi surface. In contrast to the back scattering term, the 
scaling dimension of the forward scattering term does not de¬ 
pend on KcS, which is the only Luttinger parameter modified 
by the electron-electron interaction, so the forward scattering 
term produces only local extrema (peaks) in the RKKY inter- 
actioni2i212^ Since the nuclear spin order is determined by the 
global extrema of the RKKY interaction, in what follows we 
may neglect 3 ), focus on the back scattering term, 

S'^^(rj ). Each term of Eq. (l22l) corresponds to a scattering 
process, as illustrated in Eig. [T] 



FIG. 1. Back scattering processes on the Fermi surface. indicates 
the two Dirac points with the valley index 7 = ±. The red (green) 
arrows correspond to the intravalley back scattering Ae j At,~t,rT' 
(intervalley back scattering At, ^ processes. The arrows are 

mutually independent and the spins are not shown. 
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Taking the continuum limit, we obtain 


(^) = (^) + S^Mtev,c, (»’) > (24) 

'S'&antra.a(»’) = ^l-y ’(25) 


G,a' 




'^b,inter,aM = ^ ,air)'<Pl^,a'i^) 

(7,(7' 7,^ 


.(26) 


From Eq. (1251) we see that the intravalley back scattering term, 
S'^intra Opposite for the two sublattices. We will see be¬ 
low that this term gives rise to the q = 2kp RKKY peak in 
the spin susceptibility. The intervalley back scattering terms 
in on the other hand, do not depend on the sub¬ 

lattice index and give rise to g = 2{kF — jky) peaks. The 
spin susceptibility in the absence of the feedback thus con¬ 
tains two parts, the sublattice-dependent q = 2 kF intravalley 
back scattering, and q = 2 {kF — jky) intervalley back scatter¬ 
ing, which is independent of sublattice, as in Refs. [3l53l M 
However, the latter will be suppressed when taking into ac¬ 
count the feedback (Overhauser field due to ordered nuclear 
spins), as will be discussed in Sec. HV] As a result, in spite 
of its presence in the spin susceptibility, the intervalley back 
scattering will not influence the nuclear spin order established 
by the intravalley back scattering. For clarity, we list the scat¬ 
tering processes in S^, their operators, and the locations and 
types of the corresponding extrema in the RKKY interaction 
in Table nil 

I 


Scattering process 

Operator 

Location 

Type 

Intravalley, forward 



local 

Intervalley, forward 


q ~ 2fc„ 

local 

Intravalley, back 

^1.7,aV'L7F' 

q 2kF 

global 

Intervalley, back 

^l,7,a'^L7F' 

q ~ 2(fc„ ± kp) 

local 


TABLE II. The scattering processes in S^, their corresponding op¬ 
erators, and the locations and types of the corresponding extrema in 
the RKKY interaction in the absence of the feedback. 

From now on we shall proceed with the intravalley back 
scattering term, come back to the inter¬ 

valley back scattering term when discussing the feedback 
in Sec. IIV Al Now the spin operator is expressed in terms 
of the operators '!/'^,-y,cr(r), which can be bosonized through 
Eqs. (fT^ . (fTSl l. and (fT9l l. To this end, we define the spin den¬ 
sity wave operators^^^ 

OsDW,y(.^) = XI <.(»')(27) 

(7,(7' 


such that 


OsDW,'y{''') + {pSDW,'y{''^)) ; 

(28) 

which can be written in terms of the bosonic operators, O^p 
and (pvp. 


'^Li„tra.a(^) = T X 


[‘^SI>W ,7 (’’)] ''' 
^SDW,'ri''') 


_}_g2ifcFr i g-i[^cs{r)+'r4>cA{r)+0ss{r)+'jesA{r)] _|_ ^-i[(f>cs{r)+'r4>cAir)-8ss{r)-'Y9sAir)]'\ 

27ra I J ’ 

2ikFr f 

27ra I 


g-i[if>cs(r)+'y4>cA(r)+ffes(r)+yffeA(r)] _ ^-i[rf>cs(r)+'y4>cA(r)-9ss(r)-'f8eA(r)] 


_J_ 2ikFr f 
27ra I 


g-i[^cs(r)+7't'cA(r)+^as(r)+7^aA(r)] _ ^-i[</ics(r)+7't'cA(r)-!f>as(r)-7^BA(r)] 


}. 


(29) 

(30) 

(31) 


where the Klein factors t/«, 7 ,o. are omitted because they have 
no influence. 


D. Spin susceptibility and RKKY interaction 


With Eqs. (l28l l- (l3Tl i. the spin susceptibility can be ex¬ 
pressed in terms of the bosonic fields, and calculated within 
the Luttinger liquid formalism. First, let 


= -i (S^^intrn,ai^^i)Sb,i 


intra,/9 


( 0 ) 


(32) 


where (• • •) denotes the time-ordered average corresponding 
to the electron Hamiltonian, Eq. (fTTl i. and the time argument 
appears due to the interaction representation adopted for the 
operators. In the continuum limit we have^^ 


= -2* / dte 0(f) Im 






■ (33) 


Since Eq. (O is a free bosonic Hamiltonian, the calcula¬ 
tion of the correlation functions is rather straightforward;^ 
Upon the Fourier transform, Xapi^) ~ I dre~'^'^^ 
the static spin susceptibility in momentum space reads 

Xaa(9) = -X%{<l) 

^ sin(7rg^) / 

(47r)2fiUF \27ra^ 


X 


r(i-g^)r(f-*|j(g- 2 FfcF)) 


\^^-i^{q-2KkF)) 


(34) 


where we have defined the thermal length At = 

r(a;) is the Gamma function. The spin susceptibility strongly 
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depends on the exponents, 

9oc = 9y ~ ~7 (^cS + KcA + 1“ 'Y? ) ’ (35) 

4 V ^sS KsaJ 

9 z = {KcS + Kca + Kss + Ksa) ■ (36) 

For the systems with the spin rotational symmetry, we have 
Kss = KsA = 1, which leads to 9 x = 9y = 9z and thus 
isotropic spin susceptibility, as expected. 

In Eq. (l34l i we obtained the opposite sign for x^b (?) 
plicitly. The antiferromagnetic correlation between spins on 
different sublattice sites provides a consistent picture with 
the noninteracting calculationi^i^— Our results thus consis¬ 
tently reconcile Refs. l37ll38l and [34ll53l - I?^ which obtained the 
interaction-induced boost for the transition temperature, and 
the locally antiferromagnetic coupling, separately. 

Since the RKKY coupling is related to the spin susceptibil¬ 
ity by 

(37) 

we can use Eq. (l34l i to evaluate the RKKY coupling. Its 
momentum dependence is plotted in Eig. |2] We remind that 
Eq. (l34l i contains only the contribution from the q ^ 2kp in¬ 
travalley back scattering terms, leading to the global extrema. 
The contributions from other scattering processes only give 
local extrema, and will be suppressed in the presence of the 
feedback (see Table Hill. The peak value of the RKKY interac¬ 
tion is given by 

_ A§asin(7rg^) / At 

Z2'it'^hvF \2Tra J 

which depends on the temperature through the thermal length. 
At- 


III. ANTIFERROMAGNETIC NUCLEAR SPIN HELIX 
A. Antiferromagnetic helix and magnon spectrum 

We now perform the spin-wave analysis to find spectrum 
of the low-energy excitations of the RKKY Hamiltonian in 
Eq. (|4l).i2i24 Since the RKKY interaction in CNTs is sub¬ 
lattice dependent, it leads to a different nuclear spin order 
from the ferromagnetic helical order in GaAs nanowires. We 
first consider only the long-wavelength magnons propagat¬ 
ing along the tube axis, and will include short-wavelength 
magnon excitations when estimating the transition tempera¬ 
ture in Sec. IIIIBI 

We begin by assuming that in a given cross section (i.e. 
along the transverse direction) the nuclear spins on the same 
sublattice sites in the ground state are parallel to each other 
and the spins on different sublattice sites point to the opposite 
direction, Ia(g ) = ^ helical order means these 


JabW 



FIG. 2. RKKY interaction from Eq. 1371 in momentum space. The 
interaction has peaks atg = ±2fcF with the width ~ P®' 

rameters used here are K^s = 0.2, Kca = = KbA = 1,— 

I = i, Ao = 6.0 vf = 8.0 x 10® m/s, 

kp = 4.0 X 10® m"\ a = 2.46 L = 1.0 pm, and N± = UM 
For the purpose of illustration, we choose an unrealistically short 
thermal length, corresponding to an unrealistically high temperature 
T = 100 K, to demonstrate the RKKY peaks. For realistic tempera¬ 
tures, the peaks will be much sharper. 


spins rotate within a fixed plane as one moves along the tube, 
with a spatial period TT/kp. We denote this (helical) plane as 
xy and the axis perpendicular to it as z. The confinement of 
the nuclear spins to the xy plane will be justified in Sec. El 
where we will show the modified RKKY interaction due to the 
feedback to be anisotropic; \Jap{q)\ = |Aa/j(g)| > 

We adopt the standard helical ansatz, generalized for the 
antiferromagnetic correlation between the two sublattices,— 

= aNj_I [cos{2kprj)x + sm{2kprj)y] . (39) 

As demonstrated below, this order forms the Neel order in 
a rotated basis. Although in a conventional antiferromagnet 
(i.e. Heisenberg antiferromagnet) the Neel order is not the true 
ground state, it provides a consistent basis for the spin-wave 
analysisi^^ We will make sure this is a legitimate choice here 
by checking the stability of the magnon Hamiltonian. With 
the order in Eq. ( l39b . the nuclear spins are antiferromagneti- 
cally aligned on the atomic length scale, whereas they slowly 
rotate around the helical axis (z direction) on the length scale 
of it/ kp, as sketched in Eig. [3 

To derive the magnon spectrum of the antiferromagnetic he¬ 
lix, we rotate the spin axes^— such that in the local basis 
(ej, e^, Bj) the nuclear spin model is mapped onto the Heisen¬ 
berg antiferromagnet. We write 

Ia(G) = + il{rj)y + i^{rj)z 

= ii (g )e] + il (g )ef + 7a (g )ey > (40) 

where the unit vectors in the original and new basis are related 


r(i-g^)r(f) 




(38) 
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FIG. 3. A sketch of the antiferromagnetic nuclear spin helix in the 
(a) original (x,y,z) and (b) rotated (ej,ej,ej) coordinates. The 
black and red arrows indicate the spins on the sublattice sites A and 
B, respectively. For simplicity, we do not plot the actual honeycomb 
lattice here. 


by a rotation around the helical axis 


/ 

cos{2kFrj) 

— sin(2/cFo) 

o\ 

(i\ 


sin(2/cFo) 

cos(2fcFO) 


y • (41) 

V 

0 

0 

1; 

Vv 


In the rotated basis, the Neel order forms the staggered spin 
orientation along the ej direction, i.e. = aN±Iij, as 

sketched in Fig. [3 

The RKKY Hamiltonian, Eq. (HI, becomes 

■Hrkky = “ rj)I^(rj),( 42 ) 

i,j,a,l3 fiO 


where /i,z> = 1,2,3 in spin space and the nonzero compo- 
nents of (r^ — Vj) are 

JiWi -^ 3 ) = 

= Japin - rj) cos [2kFiri - Tj )], (43a) 

Jilin -rj) = -Jipin - Tj) 

= Jipin - rj) sin [2kFin - Tj)] , (43b) 

Jilin - Tj) = - Tj). (43c) 


We now introduce the Holstein-Primakoff transformation 
for the antiferromagnet^^ 


i\ir,)=NFl-ala,, (44a) 

^Airj) = («] + “f) > (44b) 

(-«] + ’ (44c) 

Isirj) = -NfI + bibj, (44d) 

irj) = ) > (44e) 

ti('-i)=v^7 ('>; - '><). («fl 


where the higher order terms in O(iv^) have been neglected. 

Using Eqs. ( l43l l and (l44l i in Eq. (l42l i and performing the 
Eourier transform, 

“9 = (45a) 

= (45b) 

^ j 

we obtain the magnon Hamiltonian in momentum space, 

^magnon — 2iV ■E 4'Lgnon(9)25(9)^' magnon iq), (46) 

Q 

where 4/niagnon(Q) — b_q^ Ug, U_g, bq^ b—q^ . The 

8-by-8 symmetric matrix 17(g) is 


27(g) = 


-4(g) B{q) 


(47) 


^ B{q) Aiq) J ’ 
where the 4-by-4 block matrices Aiq) and S(g) are defined as 


-4(g) = 


Biq) = 


^ / 13 (g) 

0 

0 

biiq) \ 

0 

/ 13 (g) 

/i2(g) 

0 

0 

/ 12 (g) 

/i3(g) 

0 

V / 12 (g) 

0 

0 

( 13 (g) ) 


( 0 -/ 12 (g) /ii(g) 0 \ 

-/ 12 (g) 0 0 /ii(g) 

/ii(g) 0 0 -/ 12 (g) 

V 0 /ii(g) -/ 12 (g) 0 j 


(48) 


(49) 


with 


hiiq) = - [JAsiq ~ ‘^^f) + JAsiq + 2/eF) + 2J^5(g)] , 

(50a) 

/ 72 (g) = I [Jlsiq - ^kF) + Jlsiq + ^kF) - 2J1 f(9 )] , 

(50b) 
(50c) 


/ 13 (g) = 2 J^B( 2 fcF) - hi{q). 


One can check that 27(g) is positive definite, which ensures 
the stability of the nuclear spin order and justifies our ansatz 
for the antiferromagnetic helixi^^ 

The excitation spectrum of Eq. (l46t is given by twice the 
positive eigenvalues of the matrix 


-4(g) B{q) \ 

-B{q) -Aiq) )■ 

Diagonalization gives two magnon bands^^ 

/7u;W = -^j2Jl^{2kF) 

X ^ 2 J%si 2 kF) - Jlsiq - ^kF) - Jleiq + ^kF). 

(52) 

tv-o^ip = J^g(2kF) [J^F{2kF) — J^siq)]’ (55) 
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FIG. 4. Magnon spectrum of the antiferromagnetic helix. The pa¬ 
rameters used here are the same as in Fig. The blue solid and red 
dashed lines correspond to the 0 Jq^\ and energy bands, respec¬ 
tively. As in Fig. 12 here we use an unrealistically high temperature 
to illustrate the dips in the spectrum. The region of the linear disper¬ 
sion is given by the width of the RKKY peaks ~ x^'kp ’ which is 
much narrower for realistic temperatures. 


which are shown in Fig. 01 

In Fig.m one can see that there are zero-energy excitations 
at (7 = 0, ±2kp. These Goldstone modes are protected by the 
symmetries in the system; namely, the rotational symmetry 
of the nuclear spins around the helical axis and the rotation 
of the helical axis itself4^iS Around these nodes, the low- 
energy magnon spectrum exhibits linear dispersions. In the 
(nonhelical) Heisenberg model, the low-energy dispersion of 
the magnons is quadratic in the ferromagnetic case, whereas 
it becomes linear in the antiferromagnetic case.^^lH Interest¬ 
ingly we hnd that for one-dimensional helical systems, both 
locally ferromagnetio^i^i and antiferromagnetic orders con¬ 
tain low-energy magnons with linear dispersions, but differ¬ 
ences in scaling coefficients, as discussed in Sec. lIIIBl 


B. Transition temperature of the nuclear spin order without 
the feedback 

We now estimate the transition temperature by considering 
the fluctuations due to the magnons, which reduce the sub¬ 
lattice magnetization. From the magnon spectrum one can 
see that the long-wavelength magnons (g Ri 0) and their two 
replicas at g Ri ±.2kF have the smallest energy, with linear 
dispersion. In an infinite system, such excitations destroy or¬ 
der at any finite temperature. Namely, whereas the original 
Mermin-Wagner theorem^ and its extension^!! for oscillatory 
exchange interaction in free electron gas do not apply to this 
system, an extension of the theorem for a more generic Hamil¬ 
tonian, including the electron-electron interaction, rules out 
any spontaneous orders at finite temperatures in the thermody¬ 
namic limiti2i However, in a finite system with length L, the 
lowest allowed momentum is given by gi = ^ and the val¬ 
ues of the excitation momenta are discrete (not continuous). 


so that an order may be established in principle. This finite- 
size-induced finite energy of long-wavelength excitations cor¬ 
responds to a gap of the zero-energy Goldstone modes. 

In addition, here the peak in the susceptibility is so sharp 
that the region where magnons can be considered long- 
wavelength (that is, having linear dispersion) is extremely nar¬ 
row. In fact, from the analysis in Ref. [s^ it follows that for 
sample sizes realistic for nanowires and nanotubes, the long- 
wavelength magnons are completely negligible and the tran¬ 
sition temperature can be obtained by considering only the 
contribution from the short-wavelength magnons. If gi > ^ 
(equivalently kpT < which in our case is satisfied for 
any realistic length), then these magnons have approximately 
a momentum-independent energy, of order | (g = 

2fc/?)|, a property which makes the transition temperature cal¬ 
culation analytically tractable. Namely, the temperature de¬ 
pendence of the magnon occupation can be computed by 


EE' 

i—1,2 q 


1 

hu)m ^ 

e‘°BT — 1 


(54) 


where the summation over i = 1,2 includes both magnon 
bands, and the prime on the summation denotes that the Gold¬ 
stone modes are excluded. Finally, the prefactor Nj_ is re¬ 
quired to reflect the N± possibilities to flip a spin within a 
cross section for a short-wavelength magnon.— 

The order parameter, defined as the q = 2kp component of 
the normalized staggered magnetization, i.e. the normalized 
sublattice magnetization, can then be expressed as 






iti',2^ e-B 




(55) 


which equals unity for the fully ordered nuclear spin state, and 
vanishes for completely disordered phase. 

The constant magnon energies, = 

2/J^^(2feF, T)/N±, lead to a generalized Bloch law^^ 

( rp \ 3 — 2 ^^ 


with a non-universal exponent (3 — 2gx) modified by the 
electron-electron interaction. We also define 


kpTo 


Np 




(57) 


where Aa = hvp/ais the bandwidth, and 


(7(5) = isin(7rg)(27r)29 ^ 


r(i-g)r(i) 

■p f 2 — 0 \ 


2 

■ (58) 


We note that even for the noninteracting limit gx 1, the 
exponent in Eq. ( |56] | is still different from the Ta law for 
Heisenberg ferromagnet s^^'^^i^^ or the law for Heisenberg 
antiferromagnets.i^^ The nanotube radius R has two effects 
on Tq: first, as mentioned in Sec. IllB1 larger R results in less 
prominent electron-electron interaction, and thus a larger ex¬ 
ponent gx- Second, larger R, which is proportional to N± that 
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enters Eq. (fSTl i. weakens the finite-size effect, so the magnon 
occupation increases, as indicated in Eq. ( l54b . As a result of 
both of these effects, CNTs with larger R are expected to have 
a lower transition temperature. In addition, we also note that 
for CNTs not purely made of say with the concen¬ 
tration of p, Eq. (l56l l is valid upon replacement Aq pAo in 
Eq. (l52ll, so that the transition temperature will be reduced as 

To OC p2/(3-2gx)^ 

Using Eq. (fSTl) we evaluate Tq ~ 1.9 pK, which is too small 
for dilution fridge experiments. However, so far we have not 
included the Overhauser field due to the nuclear spins, which 
acts back on the conduction electrons and further stabilizes the 
order. In the next section, we take into account this feedback 
and estimate how it modifies the transition temperature. 

As a self-consistency check, we examine that the energy 
scale of the RKKY interaction, Eq. (|4]i, does not excess the 
one set by the original hyperfine Hamiltonian, Eq. (Ell. The 
former energy scale is dominated by the q = 2kp peak, 
Eq. (l38l l. Thus, after Eourier transforming Eq. (|4]i into mo¬ 
mentum space through /^(g) = we obtain 


-E'Rkky 


1 

NNl 


Jlp{q = 2kF)il{2kF)iU-‘^kF) 


NJ>XB{q = 2kF)l\ 


(59) 


where we keep only the dominant q = 2kF component, and 
I^{q = ±2kF) is replaced by its maximal value, NN±I. 
The energy scale of the hyperfine Hamiltonian, on the other 
hand, can be obtained by considering all electrons are polar¬ 
ized such that their spins locally align with the nuclear spinsi^ 
This gives 


static order, which induces a static spatially oscillating Over¬ 
hauser field that acts back on the electrons. Including the anti¬ 
ferromagnetic helix in the hyperfine coupling terms, we obtain 

. (62) 

The antiferromagnetic helix with q = 2kF gives 

(lairj)'j = aN±Im 2 kF [cos (2kFrj)x + sin (2kFrj)y ]. 

(63) 


In the continuum limit we then have 

Rfb = ^ f dr Bov,air)-Sair), (64) 

a 

where the nuclear Overhauser field is defined as 

'^Ov,air) = aBov [cos {2kFr) x + sin (2/ci?r) y\ (65) 

with Bov = AQlm2kF- The summation over a eliminates 
the intervalley back scattering terms, Eq. (l26l) . so only the in¬ 
travalley back scattering terms, Eq. (l25T l. enter TLfb- 

With Eqs. (l28T l- (l3TI) . the feedback Hamiltonian can be writ¬ 
ten as 

'Hfh = X! / + 3'<^cA + SsS + I^sa) 

7 

-f cos {(j)cS + I'PcA - dsS - idsA - 4A:_Fr)], 

( 66 ) 




Aq Uel 

2 ni 


lairj) 


NAqI 


kpa 


(60) 


where Uei = {^kp/TT)\'4’±,a\‘^ = kp/iTr'^R) is the area elec¬ 
tron density, and the electron spin is included through the fac¬ 
tor of i. Here n/ = Np/ipRa) is the area nuclei density, 
introduced in Sec. Ill Al In the second line we replace I^{rj) 
by its length, N±I, and the summation gives a factor of 2N. 
Combining Eqs. (|59ll and (gOll, we obtain the self-consistency 
condition. 


J^s('? = 2feF)< Ao^, (61) 

TT 

where JAsid = “^kp) is temperature dependent. The above 
condition is fulfilled for T = Tq. 


IV. FEEDBACK EFFECTS 
A. Overhauser field from the nuclear spin order 


where we neglected the forward scattering part because it has 
no influenceThe cosine in the second term oscillates ex¬ 
cept for the commensurate case, 2fci?a = 2p x integer. The 
commensurate case corresponds to an unrealistic gate-tuning, 
so we assume the system is incommensurate and drop the sec¬ 
ond cosine termi^i^ Consequently we have the sine-Gordon 
term 

B f 

dr [cos i(j)cs + i4’cA + dsS + ids a)] , 

7 

(67) 

which is renormalization-group (RG) relevant in the inter¬ 
acting system as discussed in Sec. II V Bl Therefore, it will 
gap out the i(j)cs + l(t>cA + dgS + I^sa) modes, but leave 
i4>cS + i4>cA — dsS — ids a) modes gapless, which can still 
effectively mediate the RKKY interaction. 

Before analyzing the Overhauser field due to the antiferro¬ 
magnetic nuclear spin helix, let us come back to the inter¬ 
valley back scattering terms of the spin operator, Eq. (|2^ . 
which would have led to a sublattice-independent ferromag¬ 
netic coupling J'XFq) = JAsid) < 0, and hence a locally 
ferromagnetic helical order with q = 2{kF — jky). 


Since Aq ^ ep, the characteristic time scales of the slow 
nuclear and fast electron dynamics can be considered to be 
decoupled. Therefore, we can treat the nuclear spin order as a 


(ifm,airj)'j = aNpIm2{kF-iK) {cos [2 (/c_f - lkv)rj] x 

+ sin[2ikF - jky)rj]y} . ( 68 ) 
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The corresponding feedback Hamiltonian is 

"Hfb,inter = ^ ^ S5^inter,a (r'j ) * ^ 


B'o. 

2TTa 


^ dr [cos {(jjcs - 7 </'sA - dsS + iQca)] , 

7 

(69) 


where B'q^ = AQlm 2 (kp-'yk^)^ the oscillating terms 
are omitted again. Here the cosine terms will gap out the 
{4>cS - i4'sA - OsS + -fOcA) modes, but leave the {4>cS + 
^(psA + SsS + I^ca) modes gapless. However, while both of 
the {(pcS + l’PsA + SsS + ldcA) modes for 7 = ± may mediate 
the ferromagnetic RKKY interaction, they generate different 
extrema at g = 2{kF T ky), respectively. On the other hand, 
for the intravalley back scattering, both 7 = ± valleys pro¬ 
duce the same extremum at q = 2kp, so the absolute value of 
the magnitude of the RKKY interaction at g = 2kF will be 
twice larger than the ones at g = 2{kF — jky). 

Comparing the two scenarios, the energy gains by form¬ 
ing these two nuclear spin orders are different because of the 
different peak heights, even though both the intervalley and 
intravalley back scattering terms lead to peaks in the RKKY 
interaction. Consequently, the ground state favors the anti¬ 
ferromagnetic helix with q = 2kF to minimize the energy. 
In addition, the gapping of half of the conduction electron 
modes reduce s the conductance by a factor of 2 , as predicted 
in Refs. l37ll!^ for materials with no valley degrees of freedom, 
which may have been observed in GaAs nanowires.— 


B. Renormalized Overhauser field 






'Ar,-h,4. 



FIG. 5. The antiferromagnetic helix corresponds to the intravalley 
back scattering processes (a), gapping out half of conduction elec¬ 
trons (b). The up and down spins are marked in black and red colors, 
respectively. The chemical potential is plotted with the blue dashed 
lines. The green arrows describe the scattering processes. The dis¬ 
persions for up and down spins are slightly shifted for clarity. 


Based on the analysis in Sec. II V Al the system will organize 
the nuclear spins to maximize the Tn 2 kp component with anti¬ 
ferromagnetic helix to lower the energy. Therefore, from now 
on we drop the intervalley back scattering contribution to the 
feedback effects, and consider only the antiferromagnetic he¬ 
lix due to the intravalley back scattering terms in Eq. ( l67l i. In 
terms of the right and left moving particles, the feedback term 
describes the {L, f) O {R, J,) scattering within each valley. 


which satisfy the commutation relation a^^'^^ 

[$^(r-i), 0 ^(r 2 )] = i|sign(r 2 - ri). (73) 

In terms of the new fields, the feedback Hamiltonian, Eq. (l67b , 
is 


V^l,7.t^R.7T + (70) 

as illustrated in Fig. |5] Notice that this is a consequence of 
the choice of helicity in Eq. ( l39] i: if the other helicity is cho¬ 
sen, namely, [cos{2kFrj)x — sm{2kFrj)y], then the antifer¬ 
romagnetic helix will correspond to the {L, 4 ,) O {R, t) scat¬ 
tering, which gaps out different spin subbands of conduction 
electrons. Even though half of the conduction electrons are 
gapped by the nuclear spin order, the feedback strongly renor¬ 
malizes the other half of the electrons, leading to stronger ef¬ 
fective electron-electron interaction, as can be seen below. 

To proceed we define a new set of bosonic fields, 

^7 = ^ [± {(pcS + IfpcA) + {dsS + idsA)] , (71) 
©7 = ^ + IfpsA) ± [OcS + IGca)] , (72) 


= dr cos (2^+), 

1 

and the electronic Hamiltonian, Eq. (fTTI) . becomes 


(74) 


u 

K L 




, (75) 


where the cross terms, such as and 

(V 0 ^) (70^), have been neglected because they are 
marginal and less important than the relevant cosine 
terms Here the modified velocity and Luttinger liquid pa- 
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rameter in the presence of the feedback are given by 



For the spin isotropic systems, K^s = Ksa = 1, we have 
u = vf/K. For noninteracting systems, K = 1 and u = vf 
are recovered. For CNTs, we have K k, 0.38 and u ~ 2.6t;i7’- 
Eqs. (TT^ and (iTSl l state that the (<i)+, 0+) and ($“,©“) 
sectors are decoupled, with the former described by a sine- 
Gordon Hamiltonian, and the latter by a free bosonic one. 
To analyze the sine-Gordon Hamiltonian, we define a di¬ 
mensionless coupling constant, y{l) = Bov{l)/^a{l) with 
Aa{l) = hu/^{l) and correlation length ^(/) = aeK Then, we 
obtain the RG flow equation for y{l)!^ 

= ( 77 ) 

where I is the cutoff length scale. In the systems under con¬ 
sideration, we always have 2 — K > 0, so y{l) grows under 
the RG flow as 

y(0 = y(0)e(2-^)\ (78) 

and the cosine term is relevant. The renormalized Overhauser 
held is then 


/.x(i-if) 

B*o, = Bo. . (79) 

The RG flow will stop when ^ exceeds the system size L, 
the thermal length At, or at a scale I*, where the coupling 
constant becomes of order 1, y{l*) ~ 1, which gives 

^(r) = ae'* = a , (80) 

with Aa = Ao(/ = 0) = hu/a. The correlation length is 
determined by the smallest scale at which any of the above 


conditions is reached. 


{ = mm{L,Ars^,«(•)}, 


( 81 ) 


In CNTs, a typical system size is of order L = 1 ym. In 
addition, Xt{T = 10 mK) Fi 1.6 mm, and ^{l*) ~ 1.8 pm, 
so T < ^{l*) -C At and we obtain ^ = L = 1 pm. 

The renormalized hyperfine coupling constant, to which the 
Overhauser held is proportional, is then 

A*=Ao(^\ . (82) 

For the noninteracting systems K = 1, and the coupling is 
the bare one. For CNTs, we get |A*| « 180|Ao| Ri 1.1 meV. 
Importantly, the renormalization is stronger than the one ob¬ 
tained within the one-band description (A* Ri 22 peV)^^^ 
because of the smaller exponent K here. With the renormal¬ 
ization we still have A* <C cf ~ 0.1 eV, so the Schrieffer- 
Wolff transformation remains well defined. The gap due to 
the antiferromagnetic helix can be obtained from the RG anal¬ 
ysis^ 


A^, = A„ 


(S)"' 


(83) 


which leads to a gap of A^ R^ 0.77 meV for our parameters. 


C. Transition temperature in the presence of the feedback 

In this section, we include the feedback into the Hamilto¬ 
nian, and compute the spin susceptibility in the presence of the 
Overhauser field. The modified RKKY interaction is propor¬ 
tional to the modified static spin susceptibility, which is now 
evaluated with the modified electronic Hamiltonian, 'Hei+'Hfb, 
where Tfei and T-Lfb are given by Eqs. (fTsT l and (l74l) . respec¬ 
tively. 

The modified correlation functions (r) = —X^gir) 
are 


eAHf)=eAiif) 

—icos{2kFr) 

2(47ra)^ 


-> z/~N —icos{2kFr) 


^ |^gtA/2<I>+(f=)g-7A/2$+(0)^ ^g7V2$;(r)g-7C2<I.-(0)^| 


(84) 


XAA(r) = 


2(47ra)^ 




+ (ev^ 


^ [$+ (f )-<i>7 (?)-e+(f=)-e- (0] g- [■i'+(o)-<i>- (o)-e+(o)-e- (o)] 


>)} 


where we used the new bosonic fields, and 0^, and de¬ 


fined r = (r, t). 


(85) 
















12 


The relevant cosine term in Eq. (f?^ tends to order the $+ 
field, which will be locked into one of the minima or max¬ 
ima of the cosine, depending on the sign of Bqv oc Aq. The 
canonically conjugated 0+ field, on the other hand, will be 
disordered. Consequently, the correlation functions of the 
field will be constants, and those containing the 0+ field will 
be exponentially suppressed4Jj2^i212^ In addition, the local ex¬ 
trema of the RKKY interaction (Table |II1) are also exponen¬ 
tially suppressed in the presence of the feedback. The physical 
picture is that since the 0i]i) sector is gapped due to the 
sine-Gordon term, its contribution to the RKKY interaction is 
much less than the gapless ($“, 0“) sector. As a result, we 
may calculate the transverse spin susceptibility by simply ne¬ 
glecting the <i>+ and 0+ fields^ and following the same steps 
as in Sec. Ill PI In comparison with Eq. (l34l i. the velocities are 
replaced by u, and acquires an extra factor of ^ because 
its first term in Eq. (l84l) contains only the gapped field 
and is thus suppressed. In addition, the essential modification, 
namely the modified exponents, is 


9x=9v = ^- ( 86 ) 

As a result, the modified static spin susceptibility for x, y com¬ 
ponents is given by 


XAAi.9) = Xaa{9) 


sin(7rga;) 

32TT^hu 



2 -2g^ 


xE 


r(i- 

3a:) r 

{¥-^ki9-2^kF)) 

2 

r 


-i^{q-2KkF)^ 



(87) 


where the thermal length now becomes At = ■ 

On the other hand, the z component of the spin suscepti¬ 
bility is exponentially suppressed by the helical order gap and 
is much smaller than the transverse component. The full ex¬ 
pression for {q) is difficult to compute because it involves 
the gapped (non-free) bosons However, since the transi¬ 
tion temperature is determined by Xa/sC?) instead of X^piQ), 
the full expression for x^^ (g) is not necessary. Nevertheless, 
to understand how the RKKY interaction depends on the gap, 
we compute the RKKY peak value at zero temperatur0S24 


Xaa{9 = 2fcT) 


1 1 
AttHu 2 — 2gz 



will be given in Sec. IV Al where the pairing gap due to the 
proximity effect is taken into account. We conclude that the 
anisotropic spin susceptibility due to the ordered spins serves 
an indirect experimental signature for the nuclear spin orderi^^ 
The feedback-modified RKKY interaction is given by 
~ The value of the RKKY peak 

mainly depends on the exponent, which depends strongly on 
the parameter KcS, as demonstrated in Eig. |6] The z com¬ 
ponent of the static spin susceptibility is exponentially sup¬ 
pressed by the helical order gap, what results in a larger 
RKKY interaction in the transverse direction |J^^(g)| > 

I J^/ 3 (g)|, further stabilizing the planar magnetic order. Ac¬ 
cordingly, this easy-plane anisotropy, in contrast to Ref. ljTlI^ 
naturally justifies the ansatz of the planar nuclear spin order. 
In Eig. |7]we plot the ratio of the RKKY peak value with the 
feedback to the one without the feedback as a function of the 
Luttinger liquid parameter, KcS- The RKKY peak is strongly 
enhanced in the presence of the feedback, and the ratio in¬ 
creases with smaller KcS, corresponding to stronger interac¬ 
tion. 



FIG. 6. Exponents and as a function of the Luttinger liquid pa¬ 
rameter, KcS ■ The blue solid and red dotted lines give the modified 
exponents Qx = jjy and in the presence of the feedback, respec¬ 
tively. The green dashed curve describes the exponent gi_i without the 
feedback. The vertical black dashed line marks the value we use to 
evaluate, KcS = 0.2. The other parameters used here are the same 
as in Fig.[^ 


We evaluate the magnon spectrum using the modified 
RKKY interaction, and estimate the transition temperature, 
repeating the procedure described in Sec. IIIIBI The tempera¬ 
ture dependence of the order parameter is the same general¬ 
ized Bloch law. 


( 88 ) 

with Am being the gap due to the antiferromagnetic helix, 
defined in Sec. IIV Bl and the modified exponent is 

For CNTs, we obtain fa 0.19, and g^ ~ 0.38. More details 
about the calculation of the RKKY peak in gapped systems 


m2kAT) = l-i^jrJ : ( 90 ) 

with a modified exponent, (3 — 2gx), and a modified transition 
temperature. 


ksTo 


2N± 



C{g.) 


(91) 
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^cS 


FIG. 7. Ratio of the RKKY peak values, Rp = = 

2kF)/ Jap{<l = 2fcF)| at T = 50 mK as a function of the Lut- 
tinger liquid parameter, KcS- The other parameters are the same as 
in Fig. 12 The black dashed line marks the value KcS = 0.2. 

which gives Tq k. 57 mK as shown in Fig. In comparison 
with the absence of the feedback, the transition temperature is 
enhanced by more than four orders of magnitude. 

The feedback not only modifies the exponents and therefore 
strongly enhances the transition temperature, but also gaps out 
the ) modes. This leaves an effective Hamiltonian for 

the gapless (<!>“, 0 ”) modes, which mix the charge and spin 
sectors of the bosonic fields [see Eqs. dlB and Con¬ 

sequently this produces a density-wave order that combines 
charge and spin degrees of freedom and reconstructs the elec¬ 
tronic statesi^iiMi The combination of the charge and spin 
degrees of freedom signifies no spin-charge separation in this 
unusual Luttinger liquid, which is equivalent to introducing a 
synthetic spin-orbit interaction.— Indeed, it has been shown 
that, upon a spin-dependent gauge transformation, a helical 
magnetic order is equivalent to spin-orbit interaction com¬ 
bined with Zeeman field . This can also be seen from the 
form of the Overhauser field, a spatially oscillating field which 
combines the spin and orbital degrees of freedom. Since the 
spin-orbit interaction is crucial for non-trivial topology, we 
now consider coupling the system to a superconductor and 
discuss the realization of MFs in CNTs. 


V. RKKY INTERACTION IN THE PRESENCE OF 
SUPERCONDUCTIVITY 

A. Spin susceptibility in the presence of the pairing gap 



0 10 20 30 40 50 60 

T(mK) 



FIG. 8. Temperature dependence of the order parameter. Top: the 
order parameter in the presence of the feedback, rh 2 kp{T), from 
Eq. (iOll. Bottom: the order parameter without (red dashed) and with 
(blue solid) the feedback from Eqs. d56b and d90b . respectively, in 
the logarithmic scale. The parameters used here are the same as in 
Fig.H 

momentum. Since there are four Fermi points in metallic 
CNTs, Cooper pairs with zero momentum can be formed ei¬ 
ther between right movers at the Kj^ valley and left movers 
at the /C_ valley, denoted as the exterior branches, or be¬ 
tween left movers at the iC-|_ valley and right movers at the 
iC_ valley, denoted as the interior branches. In the presence 
of the proximity-induced superconductivity^22.i2^ the pairing 
gaps for the exterior and interior branches are in general dif¬ 
ferent.— 

Similar to Sec. [Ill we shall first consider the system in the 
absence of feedback, and include the feedback afterward. The 
Hamiltonian now consists of the interacting conduction elec¬ 
tron terms "Hei, given by Eq. ( ITtI i. and the pairing terms 


Since the RKKY interaction in metallic phases is medi¬ 
ated by conduction electrons, it is a bit surprising that even 
in gapped phases there can still be nonvanishing RKKY peaks 
which give rise to nuclear spin orders4ii^i2^ In this section we 
show that in the presence of the superconductivity the RKKY 
interaction can still form q = zL2kp peaks, even though the 
strength of the peaks are reduced by the pairing gap. 

Here we consider only BCS-type Cooper pairs with zero 


■Hs = gs 


dr 








-(r) -F h.c. 


7,CT 



7,(T 


dr cos {'^4 >cA + - Ocs - (J^Ssa) , 

(92) 


where gg is the coupling of the pairing terms, and in the 
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second line the pairing terms are expressed in terms of the 
bosonic fields. We see that the pairing terms contain the (pss 
and 6cS fields, whereas the feedback terms contain their con¬ 
jugate fields, 9sS and (j)cs [see Eq. (IbTl il. Therefore, the anti¬ 
ferromagnetic nuclear spin helix and superconductivity com¬ 
pete with each other. The repulsive electron-electron interac¬ 
tion increases the feedback but reduces the pairing gap ."^^i^^'^° 
Nonetheless, we will see in Sec. IVI Al that the pairing terms 
are still RG relevant. 

The details of the calculation on the spin susceptibility 
in the presence of the superconductivity are given in Ap¬ 
pendix 0 The zero-temperature value of the peaks of the 
static spin susceptibility reads 


XAAi^l = 2fcF) - 


1 1 
AttHvf 2 — 2g^ 



For the noninteracting limit, -A 1, we have 


XAAi<l = 2fcF) « 


1 

4TrhvF 


In 



(94) 


which recovers the logarithmic dependence of the RKKY sus¬ 
ceptibility peak in noninteracting systems without valley de¬ 
grees of freedom4i 



FIG. 9. RKKY peak value at zero temperature as a function of the 
pairing gap in the presence of the feedback. The parameters used 
here are the same as in Fig.|2] 


B. Transition temperature in the presence of the pairing gap 

The peak value of the RKKY interaction in the presence of 

the pairing gap is given by Jg^ 2 kFi^s) = = 

2kF)\- In the presence of the antiferromagnetic helix, the 
feedback essentially modifies the exponents, There¬ 

fore, following the same procedure, the enhanced RKKY in¬ 


teraction due to the feedback can be obtained as. 


J^=2kA^s) 


Al 1 

167rAa 2 — 2gx 



(95) 


where the modified exponent gx ~ 0.19 was obtained in 
Sec. lIV cl FromEq. (i95T l and Fig.|9]we can see that the RKKY 
interaction depends on the pairing gap, or, more precisely, the 
ratio of the pairing gap to the bandwidth, A^/Aa. 

Assuming the induced superconducting gap is Ag = 0.2 K, 
the peak value of the RKKY interaction is ~ 

This relatively large peak value is due to the small ratio of 
the pairing gap to the bandwidth and small exponent gx in 
CNTs. However, as can be seen in Fig. |9] the RKKY peak 
drops quickly with an increasing pairing gap. Using Eq. (fSSl) . 
fiuim —f 2/J“^2feF/-^-L, and Eq. ( i95l l, we obtain the transition 
temperature of the antiferromagnetic helix, 

(96) 


reduced by the induced pairing gap to fcsTo Ri 17 mK. We 
note that in contrast to Eq. dMT l. the transition temperature Tq 
now scales as instead of ^ 

that the transition temperature given in Eq. ( i96b is overesti¬ 
mated because the temperature dependence of 
included in Eq. ( i95l l. However, this overestimate is negligible 
for low temperature, <C Ag. 

In the previous estimation, we assume the pairing gap of 
Ag = 0.2 K, somewhat smaller than a typical pairing gap 
(~1K) of parent superconductors used for the proximity ef¬ 
fect. Such reduction is expected due to the electron-electron 
interactions,i^2iffi and will be further discussed in Sec. IVI Al 
Alternatively, one may intentionally reduce the induced gap; 
for example, inserting a graphene sheet between a CNT and 
a superconductor. On the other hand, while a smaller Ag is 
beneficial to a higher Tq, it also results in a longer localization 
length for Majorana fermions (MFs), as discussed in Sec. I VIII 
so that a trade-off between the two parameters needs to be 
considered. 


VI. TOPOLOGICAL SUPERCONDUCTIVITY 
A. Refermionization 

We now consider the possibility to realize MFs in our sys¬ 
tem. In Refs. and 1^ it was shown that MFs may survive 
even in the presence of very strong electron-electron interac¬ 
tions, if the pairing term is RG relevant (in other words, if the 
interactions do not eliminate the pairing gap). In this case, 
the interacting bosonic Hamiltonian with spin-orbit interac¬ 
tion, Zeeman field, and pairing terms can be mapped onto a 
noninteracting fermion model with a reduced pairing gap. We 
therefore first establish that this is the case here, too. 
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To this end, we write the pairing terms in the (<1?,^ 
basis introduced in Sec. IIV Bl 

7,(7 

■ J dr [cos (—$+ + + 0 + + 0 l) 


^ 0 ^) 


= 

TTO 

which scales as 


+ cos - <&_ + 0 _,_ + 0 _)] 


(97) 


(cos (t$+ ± $_ + 0 + + 0 _)) oc 


a\ 


(98) 


with K given in Eq. (l76l l. Since the pairing gaps for exterior 
and interior branches have the same scaling dimensions, we do 
not distinguish them in Eq. (l97] | to simplify the RG analysis. 
Thus, the RG flow equation for the dimensionless coupling 

ys{l) = gs{l)/Aa{l) reads^^ 


dyail) 

dl 


2 - 1 (iT+ i 

2 1 K 


ys{i)'> 


(99) 


which gives the condition for the pairing term to be RG rele¬ 
vant. 


2-V3<K<2 + V3. 


( 100 ) 


Eor our parameters, we have K Ri 0.38, so the pairing term is 
relevant. This also justifies the reduced gap used to determine 
the order of magnitudes of the RKKY peak in Sec. IV Bl 

Here we briefly summarize the operators involved in the 
procedure. The feedback (Overhauser held), which gaps 
out the (T,t) and (i?, i) particles within each valley, can 
be described as cos(2<i)+) with the scaling dimension K. 
On the other hand, the pairing terms are written in terms 
of cos (t^’+ i + 0 j) + 01 ) with the scaling dimension 
{K + l/iT)/2. While these two terms compete with each 
other, as discussed in Sec. IV Al both of them are relevant for 
the parameters of CNTs. 

We now consider distinct exterior and interior pairing 
gapsi05ji* defined as and respectively. We then 
recast the Hamiltonian into a noninteracting fermionic model 
through the refermionization procedure. To be explicit, we 
define the slowly varying fields and 

i’'y,cr(r) = + L. 


-'7,(7'- 


^—ikpr 


such that 
and they are related to 


the bosonic fields, (p^a- and O^cr, defined in Eq. (fTSli . by 


— /K -' 

y/ lira 

= /K -' 

V 2'Ka 




A4>i<T(r)+By^(r)] 


( 101 a) 

(101b) 


After the transformation, we obtain 


Tftop = ^ j 


A,, 


-ihvFT^dr + (o'ln + 0'2T2) 

-|-As, + r72()lCr2Tl - As,_?72^20-2r2 

+Azgz(T:i\4){r), ( 102 ) 


where the 16-component spinor 0 ^(r) is defined as 


(j)\r) = 


= I r 1 (I r I r I r? 

— 1 ^_i_ -t- ? -^_i_ -b 1 ^_i_ I 1 -^_i_ I ? I"_ -t-1 _ -b 1 ^ 


-t 




-t 




,L 


'+,t’ +,4-’ +4’ -,f’ -.t’ -.i’ -T’ 

-R-r.ti -^+,4-! dj-A) ■ 

(103) 


Eurther, 5fi, tr^, and are Pauli matrices acting on 
particle-hole, valley, spin, and right/left degrees of freedom, 
respectively. A^ is the gap due to the antiferromagnetic helix 
defined in Sec. IIV Bl and As,± are defined as 


^s,± = 


Ai^) ± Ai 


(104) 


To And details how evolves with the interaction would 

require the full RG analysis. Instead, we guide ourselves by 
Refs. Isgll^ and estimate that the gap is reduced to one or¬ 
der smaller than that of the parent superconductor. Since the 
typical gap of the parent superconductor used for proximity 
effect is of order kelvin^2^^i£^ Ai^^ = A^*^ = 0.2 K are taken 
in the previous sections for the purpose of estimation. Erom 
now on we shall keep to be unfixed parameters, and 

hence As__ is nonzero in general. Einally, in Eq. (I1021 i we 
also included the Zeeman term, Az, arising from a magnetic 
held perpendicular to the helical plane (along the tube). We 
do this to break the time-reversal symmetryrSl which has in 
general profound effects on MEs. Even though the magnetic 
field along the tube also induces orbital effects,i^ 2 iS we have 
checked, by exact diagonalization, that adding them does not 
lead to any new gapped (topological) regime in the parameter 
space. Hence, we do not include such effects in Eq. (11021) for 
simplicity. 

In momentum space, the bulk Hamiltonian is characterized 
by a matrix. 


77top(fc) = hvpkTs + —^13 (ctiTi -b tT2r2) + As,+?72'5icr2'ri 
-As _772(52CT2T2 + A^pacrg, (105) 

following from Eq. (11021) upon replacing —idr -A k, which al¬ 
lows us to inspect the symmetries of the HamiltonianJS^ We 
note that while Eq. (11051 ) describes a noninteracting model, it 
retains the features of Luttinger liquid through the renormal¬ 
ized gap parameters, A^ and 


B. Topological superconductivity and MFs 

The Hamiltonian is block diagonal if decomposed into two 
pieces, "Htop = + R-lop^ with 

R-lop = \ J (106) 

where the 8 -component spinor (j)\ (r) is formed by the fields 
gapped by the nuclear spin helix, and (r) is formed by the 
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other fields, explicitly, 

-^+.- 1 ’-^-4)1 

'p\iT) = (-^+,-1-1 -^+,4,’ -^+,-1’ -^-d)- 

(107) 


of the nanotube, and is defined by the inequalities, 

(A2±A,,_)Va^-A2_+ >0, 
|Az|-|a(^)| <0, 
|Az| - aW < 0 . 


(108a) 

(108b) 

(108c) 


The corresponding 8 -by -8 Hamiltonian densities (r) and 
HfgJ (r) are obtained from Eq. (1102b . 



FIG. 10. Phase diagram on the Am-^z plane. The black solid 
curves are marked as Ci^\ and whereas the black 

dashed curves is marked as The intercepts of these curves on 

the axes are also labeled. The yellow shaded region corresponds to 
the MF regime. In the blue region the bulk spectrum is gapless. 


To find the MF solutions, we solve the Schrddinger equa¬ 
tion at zero energy while imposing the self-conjugation and 
boundary conditions on the wave functions^ The details of 
the calculation are given in Appendix iBl The results of this 
procedure are summarized in the phase diagram shown in 
Fig. m We plot it in the first quadrant of the A^-A^ 
plane, where the other three quadrants can be obtained by 
inversion symmetry about the A^ and Az axes. The for¬ 
mulas for the curves, C^\ and are given in Ap¬ 
pendix |B] Both and intersect with the A^ axis at 
Am = 'sj Ag — A^ _. and intersect with the 


Az axis at \Az\ 



|As,+ -F As,-I and \Az\ = 




I As,+ — As,_ |, respectively. The blue region corre¬ 


sponds to a nontopological gapless state. The yellow shaded 
region corresponds to a regime with two MFs at one given end 


From Fig. [To] we can see that neither the distinct exterior and 
interior pairing gaps nor the Zeeman field is necessary for this 
MF regime. We emphasize that the gap parameters, A^ and 
are modified by the electron-electron interaction, and 
therefore reflect the features of Luttinger liquid in Fig. [T0| 
The MF wave functions have composite nature and display 
multiple decay length scales, resulting in oscillations in addi¬ 
tion to the exponential decay. The localization length, ^loc, is 
determined by the largest length scale of the inverses of ki,± 
and K 2 ,e/i (defined in Appendix iBb. If the system is deep in¬ 
side the MF regime, then A^ is the largest energy scale of the 
parameters, and the localization length will be determined by 
the smaller of A^*^^ and aI*^ ; namely. 


f -, if 


< 

A« 





(^2,.)-' - , if 

Ai") 

> 

A« 






(109) 

For zero magnetic field, the localization length is the inverse 
of the smaller of the pairing gaps. 

In addition to the magnetic field along the tube, we have 
also examined that adding a magnetic field perpendicular to 
the tube, i.e. Azrj^ai, does not lead to any new gapped regime 
in the parameter space and therefore does not generate topo¬ 
logical phases with single MFs, either. Utilizing the MF wave 
functions for Az = As,_ = 0, in which case the analyti¬ 
cal solutions are available, we have checked that the MF pair 
is not mixed by a perturbation 6 H by evaluating the matrix 
element (‘hMF,i|<^ 7 f|$MF, 2 )- We found it is zero for SH cor¬ 
responding to the Zeeman field perpendicular to the tube, the 
distinct pairing gaps, the Zeeman field along the tube (includ¬ 
ing the orbital effects), and an electrostatic impurity potential. 
The fact that the MF pair is not hybridized by any of these 
terms further confirms its robustness. It is interesting to note 
that the MF pair that we find is not split in spite of the lifted de¬ 
generacy in the bulk spectrum due to the broken time-reversal 
symmetry by the external magnetic field. Through the explicit 
calculation, we conclude that Eq. (11051) supports topological 
phases with multiple MFs. It is remarkable that the criterion 
for the MFs, Eq. ( IIO 8 I 1 . is fulfilled for the parameters of CNTs. 


VII. DISCUSSION 

In the absence of experimental techniques with atomic res¬ 
olution, direct detections of the locally antiferromagnetic nu¬ 
clear spin helix are even more challenging than the ferromag¬ 
netic one, since the magnetization signal s average out to 
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zero due to the sign changes of the susceptibility between dif¬ 
ferent sublattice sitesJii Indirect experimental signatures of 
the nuclear spin order, however, can be searched for below the 
transition temperature. As discussed in the literature, these in¬ 
clude: (1) the reduction of conductance by a factor of 2 due 
to the opening of the partial (2) the anisotropic 

spin susceptibility x%p {q) Xap (?) formation of 

the nuclear spin orderj^^ (3) NMR response at the frequency 
set by the RKKY exchange due to the singular RKKY peak;— 
(4) the unusual temperature dependence of the nuclear spin re¬ 
laxation rate due to the Luttinger liquid parameters modified 
by the Overhauser field;— (5) the reentrant behavior in the 
conductance as a function of gate voltage due to the nuclear 
spin induced gapjii^ (6) the dynamical nuclear polarization at 
zero external magnetic field42. 

Furthermore, experimental probes can be implemented to 
observe the distinct pairing gaps, In general, there 

should be double-gap features below the superconducting crit¬ 
ical temperatureand the gap values should be reduced by 
the electron-electron interaction's^ Similarly, the helical gap 
due to the Overhauser field, A^, can be observed below the 
transition temperature, Tq, which decreases in the presence of 
the pairing gap. Interestingly, it has been reported that the 
NMR measurement of the double-wall CNTs, consisting of 
89% enriched inner walls and natural 1.1% outer 
walls, revealed the formation of a spin gap at low tempera- 
turesjiiS In addition, since the remaining gapless modes have 
definite helicity, CNTs may thus serve as spin filters, similar 
to the proposal in Ref. 

The localization length of MFs is set by the smaller of A^'^^ 
and Ai*^ For = 0.2-2 K and Az = 0, we ob¬ 

tain floe ~ 3-30 pm, so nanotubes with length T > 3 pm 
are needed to avoid the overlap between MFs from the two 
endsJii While increasing leads to a shorter fioc, a 

larger substantially suppresses the transition temper¬ 

ature for the nuclear spin order, so there is a trade-off between 
high To and short localization length. 

Recently, a realization of MFs in armchair CNTs driven by 
external electric fields has been proposed,— where the elec¬ 
tric fields induce the helical modesi^^iS a necessary prereq¬ 
uisite for MFs. However, those electric-field-induced MFs 
require fine tuning of the chemical potential, in contrast to 
the RKKY systems in the present work. Here, since the an¬ 
tiferromagnetic nuclear spin helix, resulting from the scatter¬ 
ing between right-moving down-spin and left-moving up-spin 
electrons, always opens a gap at the Fermi surface (Fig. |5]), 
the RKKY system does not require experimentally challeng¬ 
ing fine tuning the chemical potential. In addition, with 
the RKKY mechanism it is unnecessary to apply an exter¬ 
nal magnetic field, which is detrimental to the parent super¬ 
conductor. Further, our calculation applies to any conducting 
CNTs, and therefore does not rely on a particular chirality of 
CNTsi^ In comparison with the recently proposed spin-orbit 
coupled wires,ii2J2 nanotubes also have the advantage to 
explore MFs, owing to the aforementioned self-tuning prop¬ 
erties and the availability of high-quality samplesi22i21 On the 
other hand, since a large pairing gap reduces the RKKY inter¬ 


action and therefore Tq, parent superconductors with suitable 
pairing gaps are necessary to obtain both sufficiently high Tq 
and short fioc. 

Finally, we remark that the RKKY mechanism discussed 
here should also apply to other quasi-one-dimensional bi¬ 
partite materials, such as metallic graphene nanoribbons, in 
which hyperfine interaction is nonvanishing2i and the conduc¬ 
tion electrons mediate the RKKY interactioni^ In addition to 
the isotopically enriched materials, the antiferromagnetic he¬ 
lix can in principle be realized using magnetically doped sys¬ 
tems, where carbon atoms are substituted by magnetic atoms, 
or magnetic atoms are deposited on the material. The sub¬ 
stitutional or top-adsorbed magnetic atoms provide localized 
spins associated with a single site;^ which can take the role 
of the atoms discussed in this work, so we expect that 
the RKKY interaction can induce an antiferromagnetic helix 
in such configurations. For the plaquette (center-adsorbed) or 
bridge adatoms, on the other hand, the magnetic adatoms in¬ 
teract with an equal number of different sublattice sitesjS so 
the sublattice-dependent oscillating terms with q = 2kp in the 
spin susceptibility cancel out, provided that the couplings be¬ 
tween the magnetic adatoms and the conduction electrons on 
different sublattice sites are equal. As a result, we do not ex¬ 
pect the antiferromagnetic helix to be realized in such configu¬ 
rations. However, the remaining sublattice-independent o&ci\- 
lating terms with q = 2(fc„ ± kp) can still lead to a ferromag¬ 
netic helical orderr^li where the RKKY peaks with different 
momenta result in a beating pattern, as in two-subband quan¬ 
tum wire systems4i We also note that in graphene at half fill¬ 
ing, it was found that the plaquette or bridge adatoms lead to a 
cancellation of the oscillations in the RKKY interaction^— 
and therefore no kind of helical order will be realized in this 
case. 
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Appendix A: Spin susceptibility in the presence of 
superconductivity 

In this appendix we calculate the spin susceptibility in the 
presence of the pairing gap. As given in Sec. IV Al the pairing 
terms take the form, 

Hs = — dr cos (7(/<cA + Cr(j)sS - OcS - CTjOsa) ■ 

Tra I 

7,fT 

(Al) 

Expanding the sine-Gordon term around its minimum and 
keeping only the second-order terms, we obtain^ 

~ 2'Khvp J 
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where Ag = 2y/gJK^ is the proximity-induced pairing gap. 
In the presence of distinct exterior and interior pairing gaps, 
the pairing gap is replaced with 



2 

+ 

A«' 

1 ^o| 

> 

2 


Following Refs. l87ll95l - I^ we take the approximation u^p 
vp, and find the zero-temperature correlation functions in the 
limits of |f I = + ^ 

> ^ limit, we get 


hvp 

A, 


and 


f\ < 


hvp 

A, 


In the 


where and aI*^ are the pairing gaps for the exterior and 
interior branches, respectively. However, for simplicity we 
shall set aI*^^ = aI*^ = As in this appendix. We note that 
distinct pairing gaps are considered when investigating MFs 
in Sec.l^ 

Our goal is to recompute Eq. (|6]l with 



Xaa(at) = XAAir,T) 


— cos(2fcFr) 
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Aga 

(27ra)^ 

y/r'^ + {vf\t\ + af _ 


hvF 

.-cos{2kFr) 

a 

Uk^s+K.a) 

<1 

^ (27ra)2 

y/r^ + (t^Fh" + 


hvF 


' Ks A > 


exp 


exp 


-(CcS + Css) 

-{CcS + Csa)' 


hvp 

hvp 


, (A5) 

, (A6) 


while in the limit of |f | <C we obtain 


Xaa{t.t) = X%ir,T) = 


XAAir,T) = 


— cos{2kpr) 

(2™)^ + (u_F|r| -I- a)^ 

— cos{2kFr) 


(27ra)^ 


“I 2pa: 


29^ 


+ (u_F|r| -I- ay 


(A7) 

(A8) 


where CcS,CsS,CsA are constants of order one. For large 
distance and long time, |f| ^ which corresponds to 

small momenta and low frequencies, the correlation between 
the electrons is cut off by the superconducting gap, and the 
correlation functions exhibit an exponential decay. For small 
distance and short time, |f| ^ in contrast, the correla¬ 
tion functions retain the gapless form. If the spin rotational 
symmetry is preserved, KsS = KsA = 1, then the RKKY 
interaction is isotropic as expected. 


of the static spin susceptibility. 


Xaa(9 = ‘2.kF,uj ^ 0) Ri - 


1 


1 


AttKvf 2 — 2g^ 

A \ (2 —2s^) 


- 1 


,(A9) 


which gives Eq. (HD in Sec.lVA] 


Appendix B: MF solutions 

Eourier transforming into the momentum space and Mat- 

subara frequency domain, and taking q = 2kF, + In this appendix, we first examine the bulk spectrum of 

and w —0, we obtain the zero-temperature value of the peaks and in Eq. (11061) . and then solve the Schrodinger 
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equation for the MF solutions. The bulk spectrum of stead, we express it as the roots of the following two quartic 

defined as is too complicated to reproduce here. In- equations, 

I 


0 = 




- 2 


+ A| E^^\k) ±8A,,_A,,+Az£;«(fc) 


1 2 


+ {hvFk)'^ + 2 (A^ + Ag _ -f A^_^ — A|) {hvpky 


Az + (a.,_ + ^A2_+-A^)] [Az - (a,,_ + ^A2 + - Al) 
Az + (a,,_ - ^A2_+-A^)] [Az - (a,,_ - ^A2 + - A^) 


(Bl) 


which differ by the sign of the term linear in E^^\k). For 
nonzero Ag _ Ag ^A^, each equation gives four roots, and in 
general there are eight non-degenerate energy bands from the 
block. The bulk spectrum of El^^^ is given by EE^l. (k) 
and ±£'j-^ 2(^)’ where 


Ei%ik) 

E^lik) 


{hvpkY 


A 


(e) 


{hvpk)^ + Ai*^ 


± Az, 

(B2) 

± Az- 

(B3) 


In the absence of the superconductivity and the Zeeman field 
= aI*^ = Az = 0), half of the energy bands, E^^\k), 
are gapped by the nuclear spin helix, whereas the other half, 
E^^/l ±(.k), remains gapless as discussed in Sec.|IV] At finite 

Ag'^^A^^A^, there exists a regime where the bulk spectrum 
has band touching points at fc = Eko, where 


“ hvF + 

+2^{Al, + Al_)Al^-AlAl 


This regime is given by the inequalities, 
ko > 0; 

(a^ + a2_)a2+-a^a|>o, 


(B4) 


(B5) 

(B6) 


which is marked in blue color in Fig. [TO] In this regime 
the system is a non-topological gapless superconductor, and 
therefore not of our interest. In other regimes, the system is 
fully gapped except for the following curves, where the bulk 
gap closes at fc = 0, 


cl*' 

{Az + As,-f + Al 

- A2+ = 0, 

(B7a) 

c‘r, 

(Az-A,,_2 + A^ 

- A2+ = 0, 

(B7b) 

Cf> 

1 A^l — 1 As^_|_ + As^_ 

= 0, 

(B7c) 

cP 

\Az\ — As^+ — As__ 

= 0. 

(B7d) 

These gap closing curves are marked as C^\ C^\ 

ci‘^\ and 


and plotted as black solid curves in Fig.fTOl 


Having the gap closing boundaries, we now discuss the 
topological properties of the Hamiltonian, and investigate 
the criterion for topological phases^ To this end, we con¬ 
sider a semi-infinite nanotube with an open left end, and 
solve the Schrodinger equation at zero energy with the 
boundary condition of the MF wave function being zero 
at r = 0. Since the boundary condition is imposed in 
the real space, we need to examine it in the basis, (p = 
(c^ CA.t, CA,i, cba, Cb,i)’ which is related 

to the slowly varying fields and by Eqs. dS-dlB. 

We first focus on the Hamiltonian density of the 
first block, iTtop^(r), and solve the Schrodinger equation, 

iFtop^(r)$2^ (r) = 0. We are looking for the localized states 
at the left end of the nanotube, so we use the ansatz 
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which incorporates the self-conjugate property of MFs. This 
gives the evanescent wave functions with the exponential de¬ 
cay determined by the ki,± values. 
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which have positive real parts and thus give normalizable 
wave functions when the bulk spectrum is fully gapped. These 
Ki ± values can also be obtained by setting E^A = q and 
k = iKi,± in Eq. dB. After numerically solving the matrix 
eigenvalue equation, we find that each of the ki ± values gives 
two eigenvectors of the form in Eq. ( IB8b . which result in four 
normalizable wave functions, denoted as $22 (r), $22 
$2\('r), and $ 22 (^)- 

/ n \ 

Similarly, for the second block, we find 

four zero-energy solutions for the Schrodinger equation, 
Etop {'r)^),^yr) = 0. These are simpler and can be obtained 
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explicitly, as 
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Therefore, {‘i’±n=i 2(''’)i„=i 2(^)1’ which satisfy the 
Schrodinger equation and self-conjugation, form a set of 
eight basis wave functions. Using Eqs. dUl - (HB, 
we obtain the corresponding wave functions, B = 

{‘&± n=i 2(^)5 n=i 2 (/) basis, and examine the 

boundary condition at the left end of the tube (r = 0). For a 
given set of system parameters, (A^, As__|_, As__, A^), the 
number of MFs, Amf, is given as eight (the number of the 
column vectors in B) minus the number of the linearly inde¬ 
pendent vectors in B. 

The results are shown in Fig. [TO] The gap closing curves, 
^ and correspond to the boundaries where and 
t« 2 ,e/i vanish, respectively. Moreover, the ki_± values change 
from real or pure imaginary to complex numbers (but do not 
vanish) across the black dashed curve, marked as 


: (A^ + A2 _)A2+-A^A|=0. (B12) 
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The yellow shaded region corresponds to a topological 
regime. The MF wave functions are linear superpositions of 

2 {r) and 2 (^)’ thus display multiple de¬ 

cay length scales, arising from ki ± and K 2 ,e/i^ As a result, 
the localization length, ^loc, is determined by the largest length 
scale of the inverses of ki^± and K 2 ,e/i OT equivalently, the 

smaller of and a1*^ . 
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